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Abstract 

Multi-currency FX derivatives offer a challenging playground to the 
mathematical modelling of correlations. Quotes of liquidly traded vanilla 
options on cross FX rates, e.g. EUR/JPY, can be used to extract a great 
deal of information about the complex implied correlation structure be- 
tween the corresponding main FX rates, e.g. USD/JPY and EUR/USD. 
Including all this information in a financial model means being able to 
fit simultaneously all volatility smiles, a very demanding task. In this 
paper we propose a first solution to this problem in the class of stochastic 
volatility models. We introduce a novel multi-factor stochastic volatility 
Heston-based model, which is able to reproduce consistently typical multi- 
dimensional FX vanilla markets, while retaining the (semi)-analytical 
tractability typical of affine models and relying on a reasonable number 
of parameters. A successful joint calibration to real market data is pre- 
sented together with various in- and out-of-sample calibration exercises 
to highlight the robustness of the parameters estimation. The proposed 
model is symmetric with respect to the choice of the risk-free currency, 
opening up to new approaches for coherent pricing and risk management 
of derivatives depending on multiple currencies. 



1 Introduction 

Derivatives with multiple underlying components have attracted in recent times 
an increasing amount of attention, partly in form of public scrutiny and crit- 
icisms. Their popularity extends to the different asset classes traded in the 
financial markets, ranging from basket options written on equity stocks and 
foreign exchange rates, to CDO structures on various fixed income instruments 
and complex hybrid structures combining the different underlying types (see 

*UBS-AG, UBS, Europastrassc 1, 8152 Opfikon, (Switzerland) 

tDipartimento di Matematica, Padova (Italy) and Mathcmatisches Institut, LMU Miinchen 
(Germany) 

tDipartimento di Matematica, Padova (Italy) and Departcmcnt Mathcmatiqucs ct 
Ingenierie Financiere, ESILV, Paris La Defense (France). Acknowledgements: we are grateful 
to Jan Baldeaux, Damiano Brigo, Imran Hafeez, Eckhard Platen, Wolfgang Runggaldier for 
useful comments. 



1 



1 INTRODUCTION 



2 



e.g. |Wystup, 2006] , [Clark, 2011] , |Esqm'vel et al, 2010| ). From a client per- 
spective the appeal of such structures is clear and understandable: diversifica- 
tion reduces the exposure to individual components while the whole structure 
is cheaper to buy than multiple single-underlying options |Qu, 2010"] . The cor- 
relations between the returns of the different components are the key inputs in 
the determination of this rebate. Most often, however, these correlations are 
not observable in the market; a fact that leads to material uncertainties and 
complications in the valuation and risk management of, even rather standard, 
multi-component derivatives. 

There is nevertheless a notable and mostly overlooked example, in which 
much can be said from the market about the correlation structure. This is often 
the case of derivatives with multiple underlying FX rates (see Esqufvel et al., 20T0] ). 



Rather paradoxically, here we face the opposite problem: instead of having lit- 
tle evidence on where correlations trade, the market gives many indications 
on their structure, potentially more than any standard model, like a copula 
model, can handle. This information can be extracted from the liquid vanilla 
markets and the fact that, differently from other asset classes, appropriate mul- 
tiplication/division of FX rates are still FX rates, and hence directly observ- 
able through standard FX spot transactions. To capture the market correlation 
structure we need a model that is able to price consistently vanilla options across 
different FX pairs (see e.g. [Castagna , 2010| , |Beneder"a nd Elkenbracht-Huizi ng, 2003] ) . 
At the same time, given the high dimensionality of the problem, for any possible 
practical application we seek to retain some kind of analytical tractability. In 
this paper, we present a first proposal of such a model. 

The focus of our analysis is on the simultaneous calibration of several FX 
vanilla surfaces. This is indeed an important model prerequisite for pricing 
structures with multiple underlying FX rates. For the sake of simplicity, let 
us consider as an example the FX market composed by three currencies: say 
Euro (EUR), US dollar (USD), and Japanese Yen (JPY). All three currency 
pairs EUR/USD, USD/ JPY, and EUR/ JPY are liquidly traded, both as FX 
spot transactions and vanilla FX options. Let S d ' (t) be the spot exchange 
rate at time t as the amount of domestic (d) currency for one unit of foreign 
currency (f ) . As a consequence of the standard triangular relationship between 
the FX rates, e.g. S JPV ' EVR (t) = S' JPY ' USD (t)S' USD ' EUR (i), the price of a multi- 
dimensional option on any of the two of the three pairs is also implicitly sensitive 
to the volatility of the excluded pair. This fact is easily understood in the case 
of no-skew. To fix ideas let us consider a basket option on EUR/USD and 
USD/ JPY. The correlation between the two main FX rates Peur/usd-usd/jpy 
can be derived from the volatility of the cross via the triangular relation (see 
e.g. |Clark, 201 1| p. 228) 

2 2 2 

^EUR/JPY ^EUR/USD ^USD/JPY /-.\ 

PeUR/USD-USD/JPY — n • (.*■) 

^^EUR/USD^USD/JPY 

By varying the cross volatility <7eur/jpy an d keeping the main volatilities fixed, 
one can span over any value of the correlation between -1 and 1, yielding sig- 
nificantly different values for the basket structure. 

Similarly, if the volatility skew is included, the basket structure becomes 
sensitive to all three volatility smiles. A good pricing model must be able to 
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reproduce all three vanilla markets; theoretically, a very challenging prerequisite. 

The easiest way to proceed is to calibrate the volatility smiles of the two 
main currency pairs using typical unidimensional models and try to obtain the 
third from the correlation between them. If, however, we stick to a simple math- 
ematical modelling of the correlation, say either constant or time-dependent, we 
face two problems. Firstly, a simultaneous calibration to all volatility smiles in 
the triangle is in general difficult to obtain, as a constant correlation does not 
provide enough flexibility to match the smile of the FX cross. Secondly, the 
model is not functionally symmetric with respect to which FX pairs we choose 
to be the main ones and which one the cross. In other terms, the stochastic 
process of the cross FX rate is functionally different from the processes of the 
main FX rates it is derived from. 

This latter point, although theoretical at first sight, is of key practical im- 
portance. In principle we have the freedom of choosing any of the currencies as 
our numeraire when pricing a multi-currency structure (although some choices 
are more natural than others, e.g. the payout currency which is specified in 
the contract). A model that is able to treat symmetrically all currencies must 
be therefore preferred. Two approaches in this direction have been proposed 
in |Doust, 2007| and |Carr and Verma, 2005| . The main purpose of this pa- 
per is to specify a model with this property in the class of multi-dimensional 
stochastic volatility models. To achieve our goal, we need a paradigm change 
in the model specification. We have to forgo the standard approach of leaving 
the task of matching the volatility skew of the cross FX rate to the correla- 
tion between the main exchange rates. Instead of putting the currency pairs at 
the basis of our model, we start from the observation that any exchange rate 
may be seen as a ratio between two quantities, the value of the currencies with 
respect to some universal numeraire, and include this feature in the specifica- 
tion of the model: this reflects the point of view of the Benchmark approach 
in |Heath and Platen, 2006b] and Heath and Platen, 2006a| . In this way, our 



model does not change qualitatively depending on which perspective is used. 
The benefits of this property are manifold: i) the calibration is universal, i.e., 
it is independent on the exotic product being priced (due to the symmetry of 
the numeraire currency), ii) the price of the exotic will not vary depending on 
which perspective is chosen; moreover, Hi) the model produces symmetric and 
consistent risk sensitivities with possibly positive impacts on the way multi- 
dimensional option books can be risk managed. 

The model we present in this paper is a multi-factor stochastic volatility 
model of Heston type [Heston, 1993] , which extends the standard approach in 



Garman and K ohlhagen , 1983 . The Heston dynamics leads to an affine model 
which is known to retain analytical tractability. We will provide a complete 
discussion concerning the set of risk neutral measures and the rules to change 
between them. Rather remarkably, as a consequence of the specific Heston-type 
dynamics, the model remains functionally invariant, after parameter rescaling, 
when the risk-neutral measure is changed. This is key feature that allows ob- 
taining an universal calibration with reasonable computational effort. We will 
then test the model on real market data and show how a joint calibration of the 
volatility smiles of EUR/USD/JPY triangle is possible. In- and out-of-sample 
calibration tests will be reported to comment on the robustness of the parameter 
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estimation. 

The paper is organized as follows: we present the model in Sect. 2, initially 
using the perspective given by some kind of universal numeraire. We continue 
with the basic properties of the model, such as the positivity of the instantaneous 
covariance matrix and the presence of stochastic skewness, before presenting 
the invariance of the model and transformation rule of its parameters when 
the risk neutral measure is changed in Sect. 3. The explicit formulae for the 
characteristic function and option prices are given in Sect. 4 before turning 
to the small vol-of-vol expansions of the option prices and implied volatilities 
in Sect. 5. Finally, the joint calibration to EUR/USD/JPY market volatility 
smiles is presented in Sect. 6, together with a discussion of the procedure and 
the results, including the Feller condition and moment explosion. 



2 Multifactor Heston-based exchange model 

We consider a foreign exchange market in which N currencies are traded be- 
tween each other via standard FX spot and FX vanilla option transactions. We 
start by considering the value of each of these currencies in units of a universal 
numeraire. We will see that the discussion is independent on the exact spec- 
ification of this numeraire. To fix the ideas, one can think of it as a precious 
metal, e.g. gold, or the market portfolio like in the benchmark approach of 
|Heath and Platen, 2006b] and |Heath and Platen, 2006a] . 

Let us work in the risk neutral measure defined by the universal numeraire 
and call the value at time t of one unit of the currency i in terms of 

our universal numeraire (note that S 0,l (t) can itself be thought as an exchange 
rate, between the universal numeraire and the currency i). We model each of the 
S°' l (t) via a multi-variate Heston stochastic volatility model H eston, 1993| with 



d independent Cox-Ingersoll-Ross (CIR) components |Cox et al., 1985 , V(i) £ 



l d , and an equivalent number of driving noises, Z(t) € K d . The dimension d 
can be chosen according to the specific problem and may reflect a PCA-type 
analysis. We further assume that these stochastic volatility components are 
common between the different S '*(£)• Formally, we write 

^-^ = ( r O- r ^-(a I ) T Diag(yVW)dZ(i), i=l,...,N; (2) 
dV k {t) = K k {6 k -V k (t))dt + Z k y/VkWdW k (t), k = l,...,d; (3) 



where K k ,9 k , £fc G K are standard parameters in a CIR dynamics. Diag( x /V(t)) 
denotes the diagonal matrix with the square root of the elements of the vector 
V(i) in the principal diagonal, this term is multiplied with the linear vector 
a 1 G R d (i = 1, . . . , N); as a result, the dynamics of the exchange rate is driven 
by a linear projection of the variance factor V(i) along a direction parametrized 
by a', namely the total instantaneous variance is (a J ) T Diag(V(i))a*eft. In each 
monetary area i, the money-market account accrues interest based on the de- 
terministic risk free rate r % , 

dB l {t) =r i .B i (t)<2i, i = l,...,N; (4) 

in our gold analogy r° is the precious metal lease rate. Finally, we assume an 
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orthogonal correlation structure between the stochastic drivers 



(dZ k (t)dW h (t)) = Pk S kh dt, k,h=l,...,d, 



(5) 



together with (dZ k (t)dZ h (t)) = 5 kh dt and (dW k (t)dW h (t)) = S kh dt. 
This concludes the description of our model. 

The idea behind this approach is that each exchange rate is driven by several 
independent drivers dZ k (t) (k — l,..,d), each with an independent stochastic 
variance factor V k (i) , to which dZ k (t) is partially correlated via p k . The vectors 
a* (i = 1, . . . , N) describe by how much each of the different volatilities con- 
tribute to the dynamics of S 0,1, (t) . This correlation structure is responsible for 
the appearance of non-standard effects in the model, like a stochastic skewness, 
see subsection 12.31 

All in all, we have introduced a total number of parameters equal to Np = 
Nd + 5d (Nd from the vectors a 1 and 5 for each CIR process, K k , 9 k ,£ k , p k and 
the initial value V k (0)) to describe the volatility skew of (iV 2 — N)/2 currency 
pairs. As rule of thumb, assuming that each currency pair can be approximately 
modeled by a standard one-dimensional Heston model, which is described by 
5 parameters, around 5(iV 2 — N)/2 parameters are needed to fit all volatility 
surfaces; the value of d should be chosen to produce approximately this number 
of parameters, if not less, to avoid instabilities due to overfitting. 

Let us now turn our attention to the exchange rate S l J (t) between two 
different currencies, say i and j. We set by definition S % ^(t) :— (t)/ S°' l (t). 
By straightforward calculation, we obtain for i,j = 1, ,.,N: 



= (V - r j )dt + (a 4 - a J ') T Diag(V(t))a^ + (a* - a J ) T Diag(V V(i))dZ(i). 

(6) 



Note that at this stage we are still working under the risk neutral measure 
defined by the universal numeraire. The additional drift term in Q can be 
understood as a quanto adjustment between the currency and i. 

Proposition 1. The dynamics of the exchange rate ^ satisfies the triangular 
relation, namely for all t > and I = 1, .., N: 



This symmetry property is fundamental in order to have a model that yields 
universal calibration and consistent pricing of exotic options. In the follow- 
ing subsections, we will analyze some additional properties of the model and 
familiarize with the meaning of the different parameters. 

2.1 Vectors a 1 cannot be the canonical basis 

A rather natural choice would be to set a 1 equal to the canonical basis e l (i.e. 
the i-th element of the canonical basis of WL N , e\ = Sn, i, I = 1, .., N), then (16]), 




(7) 



Proof. Straightforward application of Ito's rule. 



□ 
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for i =/= j, reads: 

^-Ai = (r< - r*)dt + V\t)dt + JVtfjdZ\t) - y/VMdZl®, (8) 

which in the 3-currency case leads to the 3-factor Heston model. The problem 
with this choice is that the covariances (and thus the correlations) between 
different pairs are forced to be positive, 

'dS^{t) dS i ' , (t)' 



S l >i(t) ' 5*.'(t) 



V\t)dt > 0. (9) 



As we need the flexibility to define negative correlation between FX rates, the 
vectors a' cannot be taken equal to the canonical basis. These additional pa- 
rameters are needed to describe a multi-dimensional FX market. 



2.2 Stochastic covariance matrix 



By construction, the model ^ has a stochastic instantaneous (N — 1) x (N — 1) 
covariance matrix. We prove here that the specification of the model always 
leads to a positive definite covariance matrix, a fundamental prerequisite for 
any well-posed multi-dimensional model. For example, we consider the case of 
three currencies and the associated 2x2 candidate covariance matrix where we 
fix the point of view of currency i: 



( 



Cov l (t) 



dS*j(t) 



SW(t) 
/ dS l ^(t) dS^jt) 



dS^jt) riSft'ftm 

5W(t) ' S**(t) / 
dS^it) 



J 



where 



dS^jt) dS^ l (t) 



a a J 



Diag(V(t)) (a 1 - a') dt. 



(10) 



(11) 



Variances, i.e., the diagonal terms in the covariance matrix, are clearly positive 
due to the positiveness of Vk(i), k — 1, .., d. To prove that the whole covariance 
matrix is positive definite, it is sufficient to introduce the (N — 1) x d matrix B 
defined as follows: 



B, 



k = l,..,d; j = 2, 



,N, 



(12) 



j,k — «fc 

where without loss of generality we can take the point of view of currency 
i = 1. With this notation the (N — 1) x (N — 1) variance/covariance matrix 
Cov 1 (t) = Cov(t) is given by 

BDiag(V(t))B T , 



(13) 



Cov(t) 

which is positive semidefinite since for all x G K d : 

x T Cov(t)x = (S T x) T Diag(V(i))B T x > 0, 

where the last inequality follows from the positive definiteness of the matrix 
Diag(V(f)). 

In conclusion the variance-covariance matrix is well defined, and, as an im- 
portant side effect, we have the usual bound for the correlations, i.e. all corre- 
lations are bounded by one in absolute value. 
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2.3 Stochastic Skew 

To shed some additional light on the meaning of the vectors a 1 we calculate the 
skewness ? for a given exchange rate S*'-* , defined as the correlation between the 
log returns and the stochastic variance: 

w = (Noise(log g'-J), Noise(Vol(log S^))) t 

A /((Noise(log5 i J))) t v'((Noise(Vol(logS i .J))))t' 

Differently from standard single factor models, multifactor Heston models pro- 
duce a stochastic skewness. In fact, by straightforward calculation we obtain: 



= (oi-4)^V k (t)p k dt i,j = l,..,N; k=l,...,d. (15) 



Combining this term with 

(dlogS i ' j ,dlogS i ' j ) t = (a 1 - a J ) T Diag(V(t)) (a 1 - a>) dt 
and (dVk,dVk) t — £%.Vk(t)dt, we obtain 

Efc=i(°k-°*)ffc^(*)p* 

(t) = ; V ' . (16) 



ELi $V k (t)yJ(a> - ai) 1 Diag(V(*)) (a* - a*) 

This quantity is stochastic due to the presence of the variance factors V k (t) , 
k = l,..,d Christoffcrs en et al., 2009] . The vectors a* are directly related to 
the amount of skewness for each of the different exchange rates. 



3 Numeraire invariance 

Up to now we have worked under the risk neutral measure defined by our (rather 
unspecified) universal numeraire. In practical pricing applications, it is more 
convenient to change the numeraire to any of the currencies included in our FX 
multi-dimensional system. Without loss of generality, let us consider the risk 
neutral measure defined by the z-th money market account B l (t) and derive the 
dynamical equations for the standard FX rate S l '^(t), its inverse 5 J ' l (t), and a 
generic cross S^' l (t) for i,j,l = 1, .., N. 

Under the assumptions of the fundamental theorem of asset pricing (cfr. e.g. 
|Bjork, 2009| , chapters 13 and 14), investing domestic money into any foreign 
currency money-market account cannot produce a risk free return different from 
?•'. In other terms, the ratio S 1 ^ (i) (t) / B % (t) must be a local martingale under 
the i-th risk free measure Q 1 , provided that such risk neutral measure exists. 
Hence, 

S l 'i(t)B3(t)\ .S'- ■'■:/;■ «•'(/.) 



B^t) J B l {t) 

_ S lJ (t)B^(t) 
B l {t) 



a ,: - a 3 ) T Diag(V(t))a i di + (a 4 - a J ) T Diag(v/V(t))dZ(t) 
(a 1 - a-'") T Diag(yV(t))dZ^(t). (17) 
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In the last line we implicitly defined the new Brownian motion vector Z^' (t) 
under the measure Q l from the constraint of having a QMocal martingale and 
by Girsanov theorem: 

dZ{t) ql = dZ(t) + Diag(v/V(t))a l dt, i = 1, .., AT. (18) 

If we denote Q° the risk neutral measure associated with the universal 
numeraire, the Radon-Nikodym derivative corresponding to the change of mea- 
sure from Q° to Q* reads 

= exp (- £ (a*) T Diag(v^))dZ( s ) - \ J* (a') T Diag(V( S ))a^ S ) , 

(19) 

where we observe that the shift is performed only on the process Z(t). 
The Q l risk neutral dynamics of the exchange rate S l, ^(t) becomes 

dS l ' ] {t) = S l ' J {t) ((r* - r ] )dt + (a 1 - a J ) T Diag(yV(^))dZ Q, (£)) , (20) 

as desired. 

Given our assumption on the correlation structure in ([5]), we can write the 
following standard factorization under Q° 

dW k (t) = Pk dZ k (t) + yfl - p 2 k dZ£(t), fc = 1, d. (21) 

where Z _L (t) is a Brownian motion independent of Z(t). Hence the measure 
change has also an impact on the variance processes, via the correlations p k ,k = 
1, ..,d, 

dw2\t) = dW k (t) + Pk (e fc ) T Diag(yV(t))a l dt. (22) 

The component of dW k (t) which is orthogonal to the spot driver dZ(t) is not 
affected by the measure change. This choice may seem arbitrary at first sight, 
but it turns out to be consistent with the market practice. More importantly, 
Del Bano Rollin, 2008 1 shows that this assumption constitutes a sufficient con- 
dition ensuring that the foreign-domestic parity is satisfied. We finally obtain 
the dynamic equations under the new measure. With an appropriate redefinition 
of the CIR parameters 

P k =Ph, 

K k 

we can recast the variance SDE in its original form 

dV k (t) = nf(9f - V k {t))dt + Zky/W)dW$(t)' (23) 

We can now show that the Girsanov theorem used in |l8| ) has been correctly 
applied. 
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Proposition 2. The exponential local martingale is a true martingale, 

that i S E®° [|£| t ] =1. 

Proof. See the Appendix. □ 

In the Appendix we derive a proof of the previous result based on the Feller 
explosion test for diffusions as in |Wong and Heyd e7"2004| and |Mijatovic and Urusov, 2011] . 



In order to get a quick intuition, from ( 19 ) we can check that the process V has 



the same behavior at the boundaries and +oo under both measures Q° and 1 
In fact the dimension of the process V and the boundary behavior at infinity 
are the same (that is the Feller condition is unchanged and the probability of an 
explosion in finite time is zero). As a consequence, the Radon-Nikodym deriva- 



tive defined through ( 19 ) is a (true) martingale, so that the Girsanov change of 
measure is allowed. 

The invariance of the functional form of the model under measure change is 
an appealing feature of our model; other specifications of the stochastic volatility 
will almost surely break this symmetry. Financially, it makes sense to enforce 
mean reversion of the variance, other than mean explosion, yielding a condition 
on > or conversely on the original model parameters Kk, £fe, pki an d a k . 

The inverse FX rate follows from Ito calculus 



= \r 3 — r l 



(a 3 - a*) T Diag(V(t))(a* - a*)]* + (a 3 - a l ) T Diag( v /V(i))dZ Ql (t), 



which includes the sclf-quanto adjustment. Similarly, the SDE of a generic cross 
FX rate becomes 

ds"(t) s^(t) fs^m 

S*'(t) S*> l {t) \S*<j(t)J { ' 

= [r j -r l + (a? - a') T Diag(V(i))(a J - a ,; )]dt + (a* - a') T Diag(yV(t))rfZ Q ' (t). 

The additional drift term is the quanto adjustment as described by the current 
model choice. By applying Girsanov theorem again, this time switching to the 
Q J risk neutral measure, the term is removed while the CIR parameters become 

K f = K k +/>*&K-4)j 

together with the invariant p2 = p2 and ^ = Cfc ■ These are the fundamental 
transformation rules for the model parameters. 



4 Option pricing 

Together with the invariance of the model specification with respect to the 
numeraire choice, a second central feature of the model is the availability of a 
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(semi)-analytical solution for all vanilla option prices. The pricing formula itself 
is independent of the option underlying, once we work under the risk neutral 
measure associated with one of the currencies in the option and the parameters 



arc transformed via ( 26 ) 



Let us consider a call option C(S l ' J (t),K l 'i ,r),i,j — l,..,N,i ^ j, on a 
generic FX rate S l ' j (t) = exp(x^ (t)) with strike , maturity T (r = T - t 
is the time to maturity) and face equal to one unit of the foreign currency. We 
write for the CIR parameters Kk — kJ? , Ok = and so on, implicitly assuming 



that they have been transformed via ( 26 ) in the i-th risk neutral measure 



Being an affine model, the (generalized) characteristic function conditioned on 
the initial values 

ft* (w, t, t, x, V) = Ef [e w ' J (T) \x l > j {t)=x, V(t) = V] (27) 

can be derived analytically (here i = \J — 1). Standard numerical integration 
methods can then be used to invert the Fourier transform to obtain the proba- 
bility density at T or the vanilla price via integration against the payoff, with 
overall little computational effort. By applying standard arguments (see e.g. 
|Lewis, 2000| , |Sepp, 2003] ) the value of a call option can be expressed in terms 
of the integral of the product of the Fourier transform of the payoff and the 
generalized characteristic function of the log-asset pric^J 

C(S i > j (t),K i >i,r) = e- r * T ±- [ <^'(-iA,t,T,z, V)<J>(A)dA, (28) 

where 

#(A) = J e iXx (e x -K i ' j ) + dx 
is the Fourier transform of the payoff function and Z denotes the strip of reg- 



ularity of the payoff, that is the admissible domain where the integral in ( 28 ) 
is well defined. In other words, the pricing problem is essentially solved once 
the (conditional) characteristic function of the log-exchange rate is known. We 
recall the relationship between the characteristic function and the moment gen- 
crating function. In what follows we will derive the moment generating func- 
tion G l,J (w, t, t, x, V) (Laplace tranform) from which the characteristic func- 
tion is easily derived via a rotation in the complex plane <j) l ' J (cj, t, r, x, V) = 
G 1,3 {lu), t, r, x, V). It is therefore sufficient to determine the moment generat- 
ing function, as in the following proposition: 

Proposition 3. In the Multi-Heston model the conditional Laplace transform 
of the log-exchange rate is given by: 



G lJ (oj,t,T, x, V) = exp 



a 



fc=i 



(29) 



1 Hcrc wc adopt the pricing method of Lewis, 2000 who uses the characteristic function 
computed with a complex argument, also called generalized characteristic function. The 
complex argument ui typically belongs to a strip of regularity for the function (f) 1 J in order to be 
able to integrate the payoff function. On the other hand, this method avoids the introduction 
of the damping integrating factor required by the methodology of Carr and Madan ^ 1999| . 
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where for k — 1, .., d: 



A i,j M _ 2Kfc6> fc A fc ~ A fc . , .,, , , 

k [ ' ~ ° g x+ X-(r) \- A+(r)' [6U) 



2 1-e" 



B l k 3 (t) = y '- al - a{ ) — "= ; (31) 

A k = (-/s fc + ui (a{ - <) p fc &) 2 - ^ ( w 2 _ w j _ jQ 2 . (32) 

(-Kfe + uj[a k - ai) p k £ k ) ± 
A fe = " " ^ " • (33) 

Proof. See Appendix. □ 

In summary the call price is known once the characteristic function G 1 ' 3 is 
known explicitly, as in the present framework. 

5 Expansions 

As we will see in the sequel, the calibration of our model can be performed by 
relying on a standard non-linear least squares procedure, which will be employed 
to minimize the distance between the model implied volatilities and market ones. 
Model implied volatilities are extracted from the prices produced by the FFT 
routine. This procedure is quite demanding from a numerical point of view. 
An alternative approach consists in fitting implied volatilities via a simpler 
function, for example by looking at a possible relationship between the prices 



produced by the model, see (28), and the standard Black-Scholes formula for a 
suitable volatility. The next result states that it is possible to approximate the 
prices of options under the multi-Heston model, via a suitable expansion of the 
standard Black-Scholes formula and its derivatives. The proof, which is reported 
in the appendix, relies on arguments which may be found in |Lewis, 2000| and 
|Da Fonseca and Grasselli, 2011| (we drop all currency indices, it is intended 
that we are considering the FX pair). Define r = T — t and let us define 
the real deterministic functions B k °\ B k , B k , B^ k \ k = 1, ..,d as 

Bt\r)=U\-ai) J — ^ ; (34) 

\ / K, k 

B { k 1] {r) = (a* - a{) p k £ k ^ - —J ; (35) 

«£V> = (<-<) t>{—„ < 36 » 

«Pm-(4-4) foii^k — V-^r) (37) 
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and Af , A [ k l) ,A k 2) ,A { k \k = l,..,das 



2 / p-KfcT _ 1 



4 0) w = (4-4) 

4 1) w = (4-4) 3 ^a 



r e KfeT — 1 re KfcT 
— + 2 2 + 

f^k f^fo ^k 



r 2 e -K k T 

2n k 



-KkT e ~ KfcT — 1 
+ 

K k K k 



/ . \ 4 / q- P ^k * Tf> k < 

Af{r) = (al ai) 9 h ft& £ + 3^-^ + 2™ + 



2n k 



Finally, define the integrated variance as: 



-^-E(4 0, m+4 0) w^). 



(38) 
(39) 

(40) 
)■ 

(41) 
(42) 



k=l 



Proposition 4. Assume that all vol- of-vol parameters £k,k = l,..,d have been 
scaled by the same factor a > 0. Then the call price C(S(t),K, t) in the Multi- 
factor Heston-based exchange model can be approximated in terms of the scale 
factor a by differentiating the Black Scholes formula C BS (S(t), K, cr, t) with re- 
spect to the log exchange rate x(t) = lnS'(i) and the integrated variance v — <t 2 t: 

C(S(t),K,r)^C BS (S(t),K,a,r) 

d 

+ «E +B k 1 \r)V k ) d 2 xv C BS (S(t),K,a,r) 

fc=i 

+ « 2 E (4 2) (r) + B< 2) (r)V k ) d 2 vv C BS (S(t),K, a, r) 

k=l 

+ « 2 £ (4 3) M + Bf ] (r)V k ) d 3 xxv C BS (S(t),K, a, r) 



k=l 



+ 



Proof. See Appendix. 



Lfe=l 



d$ xvv C BS (S(t),K,a,r). 



(43) 
□ 



We can now state another formula, which does not involve the computa- 
tion of option prices and constitutes an approximation of the implied volatility 
surface for a short time to maturity. This formula may constitute a useful al- 
ternative in order to get a quicker calibration for short maturities. The proof is 
again provided in detail in the Appendix. 
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Proposition 5. For a short time to maturity the implied volatility expansion in 
terms of the vol-of-vol scale factor a in the multifactor Heston-based exchange 
model is given by: 



+ a 



2 ™f 



12 



fc=l / 

X: (i + 2 P i) e k (4 - 4) 4 n - ^ ( E (4- - 4) J % 



fc=l u \fe=l 



where (Tq = (a 4 — a J ) T Diag(V) (a 2 — a- 3 ) and mj = log ^^—^175 — — ) denotes 
the forward log-moneyness. 

Proof. See Appendix. □ 

6 Simultaneous calibration of the USD/EUR/JPY 
triangle 

6.1 Setup 

In this section we show an example of simultaneous calibration to three mar- 
ket volatility surfaces of options: we consider the implied volatility surfaces for 
USD/EUR, USD/JPY and EUR/JPY as observed in the FX market on a typi- 
cal day (data from 23rd July 2010), that is N = 3 with i = USD; EUR; JPY. 
The volatility sample includes expiry dates ranging from 3 days to 5 years. The 
quotes follow the standard Delta quoting conversion in the FX option market, 
we have quotes on DN, 25 Delta, 15 Delta, and 10 Delt£^] 

We try to fit simultaneously the three volatility surfaces using two stochastic 
drivers, d = 2. This choice yields a total number of parameters Np = 16, compa- 
rable to the number of parameters in 3 independent Heston models (15 parame- 
ters). This choice should not lead to overfitting instabilities. We work under the 
USD risk neutral measure to derive the option prices of the pairs EUR/USD and 



USD/JPY and the EUR measure for the EUR/JPY options, using ((28}. We cal- 
ibrate the CIR parameters in the USD measure «£ s D , 0£ SD , ££ SD , pf D , k = 1, 2. 
The parameters for the EUR/JPY are transformed to the EUR measure through 



Eqs. ( 26 ) and the invariance property of correlation and vol-of-vol parameters. 



The calibration is done via a standard non-linear least-squares optimizer 
that minimizes the total calibration error in terms of the difference between 



2 It is important to stress that in the forex market implied volatilities surfaces are expressed 
in terms of maturity and Delta (see e.g. [Wystup and Reiswic h, 2010| , |CIark, 2011| ): the 
market practice is to quote volatilities for strangles and risk reversals which can then be 
employed to reconstruct a whole surface of implied volatilities via an interpolation method 
(see e.g. |Wystup and Reiswich, 2010| , [Wystup, 20061 , |Clark, 201l1 ). Once we have the 
quotes in terms of Delta, to perform the calibration we have to convert Deltas into strike 
prices. The procedure can be found e.g. in Benedcr and Elkcnbracht-Huizing, 2003 . 
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calibrated and target implied volatities 

ElT = (Vmarkct ~ Vmodcl) ■ ( 44 ) 

n 

The use of a norm in price should be avoided as the numerical range for option 
prices may be large, thus introducing a bias in the optimization (for a more 
detailed discussion, see [Da Fonseca and Grassclli ' 201 1| ). 

6.2 Calibration results 

In Table ([l} we report the result of the calibration of the model. We performed 
the optimization considering different sets of expiries. The expiries considered 
in the largest sample are the following: 1, 2, 3, 6, 9 months and 1 year. The 
result for this particular choice of expiries is reported on the first column on 
the left. Then we repeated the experiment by excluding the largest expiry, 1 
year. The result is reported in the second column. We proceed in this way by 
excluding more and more expiries. The smallest sample is reported in the last 
column and considers only options expyring in 1 and 2 months. 

In Figs. [T] [2] and [3] we plot the market implied volatilities against those pro- 
duced by the model. The plots refer to the largest sample in Table [T] Market 
volatilities are denoted by crosses, model volatilities are denoted by circles. The 
model yields a satisfactory fit of the market points. We are not aware of any 
other stochastic volatility model that can achieve a simultaneous calibration of 
a three-currency triangle with this accuracy. 

The plots for the calibration on the sub-samples are completely analogous. 

6.3 Parameters stability tests 

In this subsection we comment on the stability of the parameters via two differ- 
ent types of analysis. We first measure the impact on the parameters resulting 
from the calibration procedure. Secondly, we fit the model parameters to a 
certain sample and then use these parameters to price an option which is not 
included in the sample. If the out-of-sample prices are close to the market, the 
model gives a reasonable description of the joint underlying FX rates dynamics. 
Moreover, the calibration can be done on a limited set of expiries, reducing the 
computation effort of the optimizer. 

As far as the first analysis is concerned, we show in Table [2] the relative vari- 
ations computed with respect to the largest sample. With the exception of ki 
we can see that there is a good degree of stability of the parameters across the 
sub-samples. Since K\ is the only value which seems to fluctuate significantly, 
we perform also a second calibration experiment, where we fix = 1, k = 1, 2. 
The results of this experiment are outlined in Table [3j The relative variation of 
the parameters can be found in Table |4j We notice that with this choice we get a 
good degree of stability, the most relevant fluctuation is now around 20% for d¥\ 



3 We do not report, for the sake of brevity, the volatility surfaces arising from this last 
experiment, but the quality of the fit is the same as before. 
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Let us now turn our attention to the out-of-sample exercise. In Tables |6j 
[7J [8] and [9] we show the difference between the market and the out-of-sample 
volatility for all sub-samples. The differences are always well below one volatility 
point, hence satisfactory. 

6.4 Feller condition 

It is well known, see e.g. [Feller, 1951[ |Andersen and Piterbarg, 2007] , that for 
square root processes represents an attainable state when the Feller condition 
is not satisfied, that is when 2nk0k < £|. In our modelling framework we 
have two volatility factors, hence we can perform the check for each factor. 
In Table [5] we report the quantity — 2^9^ — k = 1,2. We observe 
violations of the Feller condition, which constitutes a well-known fact in the FX 
derivative practice, shared with the standard one-dimensional Heston model, 
see 



Clark, 2011 



6.5 Moment explosions 

Another well established fact is that stochastic volatility models often suffer of 
pathological moment explosions which might often impact the stability of the 
pricing tools, see e.g. |Andersen and Piterbarg, 2007 , Kcllcr-Rcsscl, 2011] and 



Glasserman and Kim, 2011 . The model dynamics might lead to the explosion 



of moments, which become infinite in finite time. 

Qualitatively, the present model shares with the Heston model the same type 
of singularities. Let us define the explosion time of the n-th moment as: 

T*(n) = sup {<> : E [(S i ' j {t)) n ] < 00} . (45) 

In our 2-Heston model, the n-th moment m n (t) — E [(iS*' J (t)) ] is computed 
by ([3]) as follows: 

m n (t) =G i <i(n,Q,t,x,V) (46) 

=(5 ij '(0)) n e n ( r< - r ^ t exp{4 ,J ^)+-Bi J (i)^i(0)} (47) 

.exp{A^(t)+B^(t)V 2 (0)} (48) 

where the functions A 1 ^ (t), B^ J (f) are given in Proposition|3]with ui = n. Hence, 
the moments can be decomposed in a product of one-dimensional moment gen- 
crating functions. As a result, the moment explosion properties follow directly 
from the known results of the one-dimensional Heston model. 



In complete analogy with Keller-Ressel, 2011] (Section 6.1), from the Riccati 



ODE (60) for each factor k = 1, 2, we define the quantities Afc(n) and Xk(n) as 
follows: 

A fc (n) = (n k - u (4 - a(j p fc &) - (uj 2 - u) (a{ - a(j ; (49) 
Xk(n) = - n k + w (al - a 3 ^j p k &- (50) 
We can then apply the moment explosion formula and obtain the following 
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Proposition 6. In the multi-Heston model, 7*(n) is given by: 

1. %(n) = +oo, ifA k (n) > 0,Vfc = 1,2; 

2. T*(n) = inf fc { y=p ( arctan (^T 1 ) + ^W<o) } . */ A *(") < 
0, for some k — 1,2. 

Then we can easily calculate the time of moment explosion for all currency 
pairs. In Table [l4| we consider moments up to order 5. 



6.6 Implied risk-neutral distribution 

As a result of our calibration procedure, we have a set of model parameters 
allowing us to compute option prices which are in line with market data. The 
implied risk neutral probability densities of the different FX rates can be derived 
from option prices (cfr e.g. |Breeden and Litzenberger, 1978| ) 



(51) 

K=S 



The second derivative is approximated given a grid of option prices. We show 
the results in Figs. [3J [5] and [6] Notice that we plot the raw distribution resulting 
from prices, thus some numerical treatment is required in order to get a smooth 
distribution. All distributions feature asymmetry and leptokurtosis, features 
that are standard for returns in financial series. 



7 Conclusions 

We have introduced a new multi-factor stochastic volatility Heston-based model 
that can provide an accurate joint description of multiple FX vanilla options 
across different currency pairs. The emphasis in the model specification has 
been in the preservation of the specific symmetries of FX markets. Differently 
from other asset classes, appropriate multiplications/divisions and inversions 
of FX rates are still FX rates. Having a model that is functionally invariant 
with respect to these operations has immediate benefits: the calibrations are by 
construction universal, i.e., independent on the choice of the risk-free rate of the 
investor, leading to consistent prices and risk-sensitivities for multi-dimensional 
exotic options. 

The choice of the CIR dynamics for the stochastic variance is instrumental 
in achieving this symmetry. We have indeed proven that our model is invariant 
with respect to the choice of the numeraire once the model parameters are 
appropriately transformed. The model is always of afhne-type independently 
on which currency is used as risk free, leading to semi-analytical expression 
for all vanilla options between any of two currencies. This property is crucial 
when it comes to calibrating the model. In a standard global optimization 
algorithm we can consider together vanilla options in all currency pairs and 
achieve a simultaneous fit to the different volatility surfaces with reasonable 
computational effort. 

The model shares naturally several stylized facts with the Heston model. 
The Feller condition is often violated when fitting the model to FX volatility 
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surfaces, a common observation in the practice. Moreover, higher moments 
of the spot distribution explode at finite time; a property that might lead to 
complications/instabilities in standard numerical pricing routines. Finally, like 
any pure stochastic volatility model, our model cannot be expected to deliver a 
perfect calibration of the vanilla surfaces across all Deltas and tenors, especially 
in the short end. 

Having said that, the main result of the paper is a promising joint cali- 
bration of the model to the implied volatilities smile of the EUR/USD/JPY 
FX triangle. The fit remains satisfactory across the currency pairs, Deltas and 
tenors which were considered. Several in- and out-of-sample calibration studies 
in fact have proven the robustness of the calibration, especially once the mean 
reversion speed k has been fixed. Asymptotic expansions of the implied volatil- 
ity surface are also included as they shed light on the meaning of the different 
model parameters and can help speeding up the calibration procedure. 

The price to pay in order to obtain a consistent simultaneous calibration to 
all volatilities surfaces is that the instantaneous volatilities of the currency pairs 
do not have single dedicated drivers. Their dynamics is rather brought about 
by a linear combination of several hidden stochastic factors. As in any principal 
component analysis, it is not easy to assign a financial meaning to each model 
parameter. As this study has shown, this appealing feature has most likely to 
be traded away in order to capture the complex phenomenology of the present 
global and widely interconnected FX markets. 



8 Appendix A: Proofs 
8.1 Proof of Proposition (]2]) 

We will apply Theorem 2.1 from |Mijatovic and Urusov, 201 1| : using their no- 
tation we have = Vk and bk(x) = —pkX^ 2 a k , where we recall that by inde- 
pendency of the Vk processes the Radon-Nykodim derivative splits in a product 
and we can use a separability argument. 

If the Feller condition holds, that is when 2nkQk > then conditions a) 



and c) in Theorem 2.1 of [Mija tovic and Urusov, 201 1| hold, as from (23 1 we 
note that the boundary behavior of Vk is the same under both Q° and Q l . 

Suppose now that the Feller condition is violated. Condition a) still holds 
while condition c) does not, but the boundary I = is good. To see this let us 
compute the functions Pk{x), Sk(x) (given by formulae (16), (18) in |Mijatovic and Urusov, 201 1| ) 
in our framework: 

, 2/vfo X ^ 



p k {x) = ax 5 * exp(— 2-), 
4fc 



where a denotes a positive constant, and 



x 2 Kk e k 



Sk(x)= I p(y)dy = a I y ^ exp(-^)dy , 



where c denotes another arbitrary constant. We can now confirm using (26) in 
|Mijatovic and Urusov, 20fT] , that is good: in fact we have Sfc(0) > — oo (by 
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the assumption 2n k 9 k < £&) and the function 

(s k (x) - flfc(0))&£(a;) , , . ^ ,2 Kfe x. 

p fc (x)cr^(a;) Cfc 



is locally integrable around since s k (x) behaves like x 5 *= . Hence (26) 
holds and it follows that the Radon- Nykodim derivative ( 19 ) is a true martingale. 

8.2 Proof of Proposition (|3 



Recall that cj> (w, t, r, x, V) = G(iu),t,T,x,V), with x^(t) = log5 ij '(t). The 
functions 0, G represent resp. the characteristic function and the moment gen- 
erating function of the log-exchange rate. In order to determine these quantities, 
we first need to write the PDE satisfied by G. First of all we write down the 
dynamics of x — x 1 ^: 



dx(t) = y{r l - r j ) - i(a 4 - a J ) T Diag (V(t)) (a 1 - a J ) j dt 

+ (a 1 - a^) T Diag (-v/V(*)) dZ® (t). (52) 
We also compute the following covariation terms for k = 1, .., d: 

d(x,V k ) t = ((a' -a') T Diag(Vv) dZ<^&V^P*dZ?) 
= ((4 ~ 4) V^dZ? .fkV^PfcdZ? ) t 
= (4 - 4) V k (t)Z k p k dt. (53) 

The Laplace transform G solves the following backward Kolmogorov equa- 
tion QKaratzas and Shreve, 199"!] ): 

d 2 G i s ,\„ . 1^9 2 G 



fe=i fc=i K 

1 - r 1 ') - - (a 1 - a J ) T Diag (V) (a 1 - 



dG 
dx 

d BG 

+ J2gv K k(e k -V k )=0 (54) 
fc=i k 

with terminal condition G (cj, T, 0, a;, V) = e ux with w € M. In order to solve 
this problem we look for an exponential affine solution of the form: 



G (u, t, t, x, V) = exp ^A(t, T) + £ B fc (t, T)V k + C(t, T)x^j , 



(55) 



for some deterministic functions A,B k ,C that may depend on both t, T. Upon 
substitution of the guess and recognition of the terms we obtain the following 
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system of d + 2 ODE's: 



dA 
~dt 



Y,B k (t : T) Kk e k + (r l - ft) C(t,T) = 0; 



fc=i 



1 

dt 



+ -C 2 {t, T) (4 - 4) + C(t, T)B fc (t, T) (a{ - a{) p k £ k 
l-B 2 k (t, T)£ - I (4 - ai) 2 C(t, T) - B k (t, T)n k = 0; 



dC 
dt 



= 0, 



(56) 

(57) 
(58) 



with terminal conditions: A(T,T) = 0, B k (T, T) = 0, C(T,T) = uj for 
= l,..,d. From (58 1 and its terminal condition, we deduce that C(t,T) = lj 
for t G [0,T], so we can rewrite the system as follows: 



dA 
~dt 



J2 K k 6 k B k {t, T) + (r l - r j )w = 0; 



fc=i 



5S fc 1 



0f 



B£(t, T)^ + -« fe + w 4 - 4 P*& r ) 



(59) 



(60) 



Now for d = l,..,d we assume that B k (t,T) can be written by means of a 
function E k (t,T) and set: 



B k (t,T) 



d E k (t,T) 



tit 



then the solution for (60 1 is: 



B k (t,T) 



A+ e -v^(T-t) _ A - 



with 



A 



At 



fe = (-«* + w (a l fe - 4^ pfc^) - Cfe ~ w) (al - 4) 
-K fc + w (a£. - 4) Pfc^fe) ± V^k 



(61) 



(62) 



(63) 



(64) 



Equipped with the solution for B k (t, T) we can now compute A(t, T) as follows: 

,T Q 

du 

fT d 



A(T, T) — A(t, T) = 



A(u, T)du 

A(t,T)= [ V" K k 8 k B k (u, T) + (r* - r°) udu 

d T 

(T-t) + y^K k 6 k / B k (u,T)di, 



[r — r J u 



2/c fe fe 



(r*-r'>(T-i) + £ ,2 
fe=i 



9/ 



£fc(*,T) 



E k (t,T) 



du, (65) 
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which implies that the solution for A(t,T) is 



A(t,T)=(r < -W>(T-t) + ^^log 



2n k 9 k , — A fc 



f. \ + e Ar(T-t) \- A7"(T-t) 

fc=l A fc e fc V ^ - A fe e fc V ; 



(r'-r^(T-() + E^( f ' T ). (66) 



fe=i 



where the functions A k (t, T) are implicitly defined by the last equality for 
d = f, ..,d. Now we obtain the statement of the proposition once we replace 
B l k 3 (T) = B k (t, T), A)f(r) = A k (t, T) with r = T-t. 



8.3 Proof of Proposition (El) 



The starting point is given by the Riccati ODE (60) expressed in terms of 
time-to-maturity r = T — t and perturbed by introducing the vol-of-vol scale 
parameter a: 



dr 



^B 2 k {r)a 2 e k + (-Kfc + w (4 - Ofc) Pfcaffc) ^(r) 
^(4-4) 2 , fe = l,..,d. (67) 



We consider the following expansion in terms of a: B k (r) = B k ,o(T)+aB k ,i(T) + 
d 2 B kt 2{T). By plugging in the expansion and upon recognition of terms we 
obtain the following system of ODE's: 



2 2 

= -K k B k , (r) + (4 - 4 ) ; (68) 



dB ki0 

dr 
dB kA 

dr 

dBk > 2 - _, . n. .,-> , /O _ „' 



-Kfc-B fcj i(r) +w (<4 - ajQ p k £kB kfl {T); (69) 
(Vr -«fcS fe>2 (r) + w (4 - 4) ft&B M (r) + ^, (r)^« (70) 

2 _ 

If we denote 7 := ^ 2 ^ then the solutions are easily computed as: 

B k , (r) = B^o^e-^ + e-^ r J e Kfctl 7 (a{ - <) ' 
=0 

= 7^ 0) (r); (71) 

B fc ,i(T) = B M (0) e- K ^ + e -"* T ^V^w - ojQ p fc & 7 B< 0) (u)du 
=0 

= "7£i 1} (T); (72) 

= B M (0) e- K » T + e~ K * T J e^ u u? 1 (a\ - a(j p k £ k B k 1] \u)du 
=o 

+ e- KkT / e Kfcil 7 ^ (b£ 0) (u)) du 

= c 2 7 6f(r)+7 2 ^ 2) W. (73) 
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Then we can write the function B k (r) as follows: 

B k (r) = 7 fif (r) + a^B[}\r) + c? [u?^\r) + T^fW) ■ (74) 



A direct substitution of (74) into (59) allows us to express the function A(r) 

A(t) = lj (r l - r j ) t + V K k 9 k / B k (u)du 
k=i J ° 

T + jJ2 Kk9k I B ( k 0) (u)du+u }1 aY,^kO k [ B?\u)di 
Jo t _ n . Jo 



■ W 



it . (75) 



We consider then the price in terms of Fourier transform as in ( 28 ) by replac- 
ing the argument uj — iA. A Taylor-McLaurin expansion w.r.t. a gives the 
following: 



C(S(t),K,r) 



2tt 



e iA(r'-H)r+iAx+ 7 ELl(<"M+< ) M^)$(A)dA 



fc=l 

2tt 



/iAe 1A ( r ^ r ^ T+iA2:+ ' rl: "= 1 K n)(r)+ ^ 0,(T)V ' c ) 



$(A)dA 



fc=i 

-7**t 



2tt 



7 2 e 



3 iA(r--H)r+iAx +7 i:L 1 (-4i 0) M+ei 0) M^) $(A:)dA 

a 2 E(4 3) (r)+Bi 3) (r)14) 



fc = l 



2tt 



a 

T 



Lfc=l 



2tt 



7iA 2 e iA(r i -^)r+iAx+ 7 EL 1 K ) W+4° ) WV fe ) $(A)dA 

7 2 iA 2 e iA ( r, - rJ > +iAl+7ELl K >(T)+ ^ 0)(T)yfc )$(A)dA. 



(76) 



Recall now from (42 ) the definition of the integrated Black-Scholes variance. In 
the previous formula, in the first term on the right hand side, we recognise the 
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Black-Scholes price in terms of the characteristic function when the integrated 



variance is v = o r: 



C BS (S(t),K,a,r) 



2- 



e iA(r i -^)r+iAx+<iiV^^(A) ( iA, (77) 



so that the price expansion is of the form: 
C(S(t),K,r)^C BS (S(t),K,a,r) 



+ aJ2 (4 1} W +B£ ) (r)V k ) d 2 xv C BS (S(t),K,a,r) 
fc=i 

d 

+ a 2 ]T (4 2) (r) + (r)V k ) d 2 vv C BS (S(t),K, a, r) 

k=l 
d 

+ « 2 E (4 3) to + B ( k 3) (r)V k ) d 3 xxv C BS (S(t),K, a, r) 



k=l 



Lk = l 



d xxvv C BS (S(t),K, cr,r) 



(78) 



From the previous expression we can deduce the relation defining the deter- 



ministic functions B k h \A k h \ h ~ 0, 3. 



BjV)=(»l-«0 



2 1 _ e -* k T 

B k 1] (r) = (4 - 4) Pktke-^ f\^Bf\u)du- 



B ( k 2) (r) 



2n k 



f 

n k r e KkU (B { k Q) {u)y 



du: 



Bi 



»(r) = (4 - 4) p k £ k e-^ f er*"B$\u)du 



and 



Af(T) = n k 6 k I B k u '(u)d 



(79) 
(80) 
(81) 
(82) 

(83) 
(84) 
(85) 
(86) 



8.4 Proof of Proposition (J5J) 

We follow the procedure in |Da Fonseca and G rassclli , 2011| . We suppose an 
expansion for the integrated implied variance of the form v = of mp T = Co + Q ; Ci + 



A < k 1) {T)=K k 9 k I B K k \u)du; 



A { k \r) = K k Q k I B k z> {u)d 



A k 3 \r) = n k e k I B k *>{u)du. 
Jo 

Computing the trivial integrals completes the proof. 
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a 2 (2 and we consider the Black-Scholes formula as a function of the integrated 
implied variance and the log exchange rate x = logS*: C BS (S(t), K,<j,t) — 
C B s {x(t), K, <r 2 mp T, r). A Taylor-McLaurin expansion gives us the following: 



C BS (x(t),K, a 



imp 



= C BS (x(t),K, Co, t) + aCid v C BS (x(t),K, Co, t) 

2 

+ y (2( 2 d v C BS (x(t),K, Co, r) + Ci 2 ^ 2 2 C BS (x(t), K, Co, r)) . 

(87) 



By comparing this with the price expansion ( 78 ) we deduce that the coefficients 
must be of the form: 



Co = vq; 



ELi [A^{r) + Bt>{r)V k )dl v C B ^ 



?(!)/ 



(89) 



-Ci 2 ^C BS +2Eti (4 2) (r)+B[ 2 \r)V k ) d 2 vv C BS 



n4 /o 



2<9, J C R 



(90) 



where the Black-Scholes formula C BS (x(t),K, <J 2 mp r 1 r) is evaluated at the point 
(a^-JT, Vq,t). In order to find the values of Ci,C2, we differentiate (34)-(37) thus 
obtaining the following ODE's: 



dB 



(0) 



9r 



dB 



(i) 



0r 



as 



(2) 



Or 



dB 



(3) 



5r 



- Kfc ^ 0) (r)+(^- fl i) 2 ; 

-« fc Bi 1) (r) + (4-4)^Bi 0) W; 



-^Bi 2) (r) + i^Bf(r) 2 ; 



- Kfc sf ) M+(a 2 fc -ai)p fc ^i 1) (r). 



We consider a Taylor-McLaurin expansion in terms of r: 



^ 0) (r) 



Si 2) (r) 



2 



K fc a k - a% 



4 X) W = y(4-4) ^-§"^(4-4) 



3 (4-4)' 



(91) 
(92) 
(93) 
(94) 
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Noting from ( 38 1-( 41 ) that are one order in r higher than the corresponding 

(i) 

B k , the following approximations hold: 



E (4V) + ^ 0) (r)F fe ) = E (4 - 4) V^r + o(r); 
fc=i fe=i 

E W + B« (r)7 fc ) = E P*& (at - 4) " V k y + o(r 2 ); 
fe=i fc=i 

£ (4 2) (r) + f(r)%) = E 4 2 (4 - 4)^4 + °(A 

k=l k=l 

E U 3) (r) + Bf (r)%) = E PUl (4 " 4) * ^7" + «(- 3 )- 



fe=i 



fc=i 



(95) 
(96) 
(97) 
(98) 



We introduce two variables: the log- forward moneyness my := log I g J 

and V" := (a 4 - a J ') T Diag(V) (a 1 - a 3 ') r. Then, from |Lewis, 2000] , we con- 
sider the following ratios among the derivatives of the Black-Scholes formula: 



(99) 
(100) 
(101) 
(102) 



d 2 v C BS (x,K,V,t) 


1 


my 


d v C BS (x,K,V,t) 


= 2 + 


V ' 


d 2 v C BS (x,K,V,t) 


o 
f 


1 1 


d v C BS (x, K, V, r) 


= W 2 


~ 2V ~ 8' 


dl xv C ss {x,K,V,t) 


1 


m f -l m 2 


d v C BS (x,K,V,t) 


= 4 + 


V V 2 ' 


di xvv C BS (x, K, V, r) 


4 

_ / 


m 2 (my — 1) 


d v C BS (x, K, V, t) 


2V A 


2V 3 



Upon substitution of ( 95 1-( 102 ) into (88)-(90), we obtain the values for Q,i — 
0, 1, 2 allowing us to express the expansion of the implied volatility. 

Co = (a 1 - a') T Diag(V) (a* - a') r; (103) 



my 



c 2 = 



■nr. 



(a* -a?) Diag(V) (a* - a?) 



(a* -a?) Diag(V) (a' - a-?') 

s(gs(4-«i)V k ) 



(104) 



16 ( a * - a?) Diag(V) (a 4 - aJ) 



(105) 
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Table 1: This table reports the results of the calibration of the model. We 
concentrate on the two factor case. For each column, a different number of 
expiries, ranging from 6 to 2, is chosen. More specifically, 6 means that the 
following expiries are considered: 1, 2, 3, 6, 9 months and 1 year, whereas 5 
means that the longest maturity, i.e. 1 year is excluded from the sample. We 
proceed analogously in the subsequent columns by excluding the longest expiry 
date up to the point where we perform the calibration on the 2-sample, where 
we fit the smile at 1 and 2 months. We consider market data as of July 23rd 
2010. The reference exchange rates are EURJPY=112.29, EURUSD=1.2921 
and USDJPY=86.90. Res. norm, is the residual of the objective function for 
the given set of parameters. 
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Table 2: In this table we consider the calibration on the largest sample as a 
basic case. We report the percentage difference between the model parameters 
resulting from the subsamples 
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Table 3: This table reports the results of the calibration of the model. In this 
case we are assuming Kk = l,k = 1,2. For each column, a different number 
of expiries, ranging from 6 to 2, is chosen. Res. norm, is the residual of the 
objective function for the given set of parameters. 
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Table 4: In this table we consider the calibration on the largest sample as a 
basic case, when = 1, k = 1, 2. We report the percentage difference between 
the model parameters resulting from the subsamples 

6 5 4 3 2 



k = 1 -0.1715 -0.1841 -0.2076 -0.2276 -0.2445 
k = 2 -0.6745 -0.6640 -0.5086 -0.5153 -0.5768 



Table 5: For all k = 1, 2 and for each sample we report the quantity 2n^Q^ — 
In all cases the quantity is negative and its absolute value is a measure of the 
violation of the Feller condition. 
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Table 6: This table reports the raw difference between the market implied 
volatility and the implied volatility for 1 year, when we calibrate the model 
to the previous 5 expiries. Blanks on the first column reflect missing market 
data for 15DC and 15DP. 
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Table 7: This table reports the raw difference between the market implied 
volatility and the implied volatility for 1 year and 9 months, when we calibrate 
the model to the previous 4 expiries. Blanks on the first two columns reflect 
missing market data for 15DC and 15DP. 
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Table 8: This table reports the raw difference between the market implied 
volatility and the implied volatility for 1 year, 9 and 6 months, when we calibrate 
the model to the previous 3 expiries. 
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Table 9: This table reports the raw difference between the market implied 
volatility and the implied volatility for 1 year, 9, 6 and 3 months, when we 
calibrate the model to the previous 2 expiries. 
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Table 10: This table reports the raw difference between the market implied 
volatility and the implied volatility for 1 year, when we calibrate the model to 
the previous 5 expiries and nu — 1, k = 1,2. Blanks on the first column reflect 
missing market data for 15DC and 15DP. 
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Table 11: This table reports the raw difference between the market implied 
volatility and the implied volatility for 1 year and 9 months, when we calibrate 
the model to the previous 4 expiries and Kfe = l,k = 1,2. Blanks on the first 
two columns reflect missing market data for 15DC and 15DP. 
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Table 12: This table reports the raw difference between the market implied 
volatility and the implied volatility for 1 year, 9 and 6 months, when we calibrate 
the model to the previous 3 expiries and Kk = 1, k = 1, 2. 
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Table 13: This table reports the raw difference between the market implied 
volatility and the implied volatility for 1 year, 9, 6 and 3 months, when we 
calibrate the model to the previous 3 expiries and Kk = 1, k = 1,2. 



Order 


gUSD.EUR 


g.IPY.USD 


g.JPY.EUP.. 


1 


+ 0O 


+0O 


+ 0O 


2 


+oo 


12.1962 


5.1612 


3 


+oo 


2.9537 


2.0580 


4 


3.3968 


1.7990 


1.3763 


5 


2.0070 


2.0819 


1.0614 



Table 14: Times of moment explosions for moments up to order 5 for the three 
currency pairs. 
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Figure 1: Calibration of USD/EUR implied volatility surface. 
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Figure 2: Calibration of USD/JPY implied volatility surface. 
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Figure 3: Calibration of EUR/JPY implied volatility surface. 
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Figure 4: Implied risk- neutral conditional distributions for the USD/EUR ex- 
change rate inferred from option prices. 
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Figure 5: Implied risk-neutral conditional distributions for the USD/JPY 
change rate inferred from option prices. 
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Figure 6: Implied risk- neutral conditional distributions for the EUR/JPY 
change rate inferred from option prices. 
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